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Abstract 

We discuss universality properties of blow-up of a classical (smooth) solutions 
of conservation laws in one space dimension. It is shown that the renormalized 
wave profile tends to a universal function, which is independent both of initial 
conditions and of the form of a conservation law. This property is explained 
in terms of the renormalization group theory. A solitary wave appears in log- 
arithmic coordinates of the Fourier space as a counterpart of this universality. 
Universality is demonstrated in two examples: Burgers equation and dynamics 
of ideal polytropic gas. 
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1 Introduction 

Conservation laws constitute the foundation of modern physics. They are linear when 
describing single or several elementary particle systems, and become nonlinear for 
macroscopic systems, motivating the vast research area of nonlinear wave dynamics. 
The important class of conservation laws corresponds to idealized models, when the flux 
functions depend on state variables only (neglecting derivative terms corresponding to 
dissipation etc.), see, e.g., [HID]. in these models, smooth initial conditions typically 
give rise to a blow-up (singularity) in finite time followed by formation of a shock wave. 
Development of such finite-time singularities is a classical subject described in every 
book on nonlinear waves, e.g., [TJ IH [TO] . 

It is well understood that scalar conservation laws in one space dimension capture a 
lot of what happens in general systems. In this respect, the inviscid Burgers (or Hopf) 
equation 

du du 
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plays a special role. Solution of this equation is easily constructed by the method 
of characteristics. Each characteristic carries a constant value of u, and the blow-up 
occurs when characteristics cross. 

In this paper we show that the wave profile u(t,x) creates a universal "core" just 
before the shock formation, as a consequence of a cusp catastrophe in (t, x, u) space. 
This core is described, up to scaling symmetry, by an analytic function, which is inde- 
pendent of initial conditions as well as of the flux function. We explain the universality 
using the renormalization group approach, and discuss its relation to more sophisti- 
cated universality phenomena described by the theory of renormalization group like 
critical phenomena of second-order phase transitions [IT] among many others [6j |8]. 
Finally, we show that in logarithmic coordinated of the Fourier transformed function 
u(t,k), the blow-up is mapped to a stable solitary wave moving with constant speed. 
As an example, we describe universal structure of shock formation in ideal polytropic 
gas. 



2 Universal structure of blow-up 

First, let us consider the inviscid Burgers equation ([I]) with smooth initial condition 
u(to,x) = uq{x). Its solution u(t,x) constructed by the method of characteristics has 
the implicit form 

x = x + u (x )(t -t ), u = u (x ), (2) 
where Xq is an auxiliary variable. For spatial derivative of u(t,x), we have 

du _ du/dxp _ u' (x ) 

dx dx/dx Q 1 + u' (x )(t — t ) ' 

The denominator vanishes at t = t — 1/u' (xq). This yields the well-known result that 
the classical solution blows-up along the characteristic with the minimum negative 
value of u' (xo), followed by the formation of a shock wave in a weak solution. 

We choose the origin of time and space so that the blow-up singularity appears 
at t — x — 0, and consider the classical solution in the interval to < t < 0. Also, 
we can take u = at the singularity, which can be achieved by the transformation 
x (->• x — uo(0)(t — to) and u (0), which leaves (pQ) invariant. In this case, the 

initial condition uq(x) satisfies 

wo(0) = 0, t = lK(0)<0, <(0) = 0, <'(0)>0, (4) 

which are the blow-up conditions at t = x = u = 0. Using (j3J) in (j5J), for small xq, we 
obtain 

x = ut- 6 ° Xq + o^), u = u' (0)x + o(x ). (5) 

Equivalently, 

u"'(0) 

x = ut-cu i + o(u i ), c=-^-L- 4 >0. (6) 



In particular, at t = 0, we obtain the well-known singular dependence u ~ —x 
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Figure 1: (a) Convergence of scaled profiles u\(to, x) (thin black curves) to the universal 
function w(t Q ,x) (bold red curve) for the inviscid Burgers equation (Tj[|). The initial 
condition is u (x) = — (y / 7r/2)erf x at to = — 1, and A = 1, 10, 10 2 , 10 3 . (b) Similar 
convergence results ( fl5l) for blow-up of a simple wave in ideal polytropic gas. The bold 
red curve represents the scaled universal function w(to,x)/ f"(0), and thin back curves 
are the scaled profiles U\(t ,x) = Q\U(t ,x) corresponding to the density variation U 
determined by (jIE]) ? (jTH]), (E2T]) . 



Consider now the renormalization function 

Ux (t, x) = g x u{t, x) = A 1/3 u( A~ 2/3 £, X^x). (7) 

It is easy to see that u\ is a new (scaled) solution of (CQ), so (CQ) represents one of the 
symmetries of (CQ). Multiplying both sides of (jEj) by A and making the substitution 
t (->• A~ 2//3 t, x I—)- A _1 a; yields the equation for u\(t,x) as 



x 



MA t- CM 3 + Ao((A~ 1 / 3 MA ) 3 ). 



The last (correction) term contains powers (\^ 3 u x ) n with n > 3. Hence, it vanishes 
in the limit of large A, and (J5]) takes the exact form 

x = wt — cw 3 (9) 

for the limiting function 

w(t,x) = \im u\(t, x). (10) 

Equation ([9]) determines the universal function w(x,t) for t < 0, which is indepen- 
dent of initial conditions. Note that different values of the coefficient c correspond to 
w(t,x) determined up to scaling y/cw(t,x/y/c), which is a symmetry of the Burgers 
equation (TT]). The function w(t,x) is a self-similar solution, 

Q\w = w. (11) 

A numerical example of convergence of u(t, x) to the universal function w(t, x) is shown 
in Fig. QJi. 

The same function w(t,x) describes the limiting universal behavior for a general 
scalar conservation law 

dU df(U) , , 
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with an analytic flux function f(U) and initial condition Uq(x) at t = t . This fact can 
be shown as follows. Solution of equation (1121) has the form 



x = x + f'(U (x ))(t - t ), U = U (x ), (13) 

which reduces to solution (El) of the Burgers equation by the substitution u = f'(U). 
We assume that the coordinates for time, space and state are chosen such that (jlj) 
holds for Uo(x) = f'(Uo(x)), and u = f'(U) is locally invertible as U — g{u) with 
g(0) = /'(0) = 0. This implies the blow-up at t = x = with U = u = 0. We will use 
the relation 

g x u n (t, x) = \ 1/3 u n (\- 2/ H, A _1 x) = A (1 ~ n)/3 <(t, x). (14) 

Expanding g(u) in Taylor series and using (1101) . (1141) . we see that all terms with n > 1 
vanish for large A. Since g(0) = and g'(0) = l//"(0), the remaining terms yield 

lim g x U{t, x) = lim g x g{u(t, x)) = w(t, x)/f"(0). (15) 

A— >oo A— >oo 

Conditions (1101) and (ITS]) with the function w(t,x) given by demonstrate strong 
universal character of the blow-up of a classical solution for a scalar ID conservation 
law. We see that, when the singular point is approached, a part of the wave profile 
u(t, x) takes an absolutely universal form independent both of the initial condition and 
of the flux function (up to the scaling transformation). In particular, derivatives of all 
orders turn out to be universal for the scaled solution near the singular point. 

Note that expressions fl6T)- f[TUl) relate the blow-up with a cusp catastrophe in the 
space (t,x, u), see [2], so that the universality of the blow-up represents a general 
property of this catastrophe. Since the cusp catastrophe is a generic phenomenon, the 
universality results are extended directly to systems of coupled conservation laws in 
one space dimension, for which a general solution similar to (TT3"j) does not exist. 

As an example, let us consider formation of a shock wave in a simple wave solution 
for one-dimensional flow of ideal polytropic gas. The density p(t, x) in this wave is 
described implicitly by 

x-Xp 7+ 1 rr- ( 7 -i)/2/ x ( s na s 

-, — t = tV a 1Po (xo), P = Po(x ), (16) 

t - 1 7-1 

where po(x) is the initial condition at t = to, see, e.g., [3]. We will use the values 
7 = 5/3, A = 3/5, and po{x) = 2 — arctanx. Then expressions (fT6l) take the form 

x = x + 4(2 — arctan^o) 1 / 3 ^ — t ), 
p = 2 — arctanxo, 

which can be written in the form (1131) for the density variation in moving frame 

U (t, x) = p(t, x + x 1 + v l {t- to)) - pi (18) 

with 

f'(V) = 4(U + p 1 )W-v 1 , 
Uq(x) = 2 — arctan(x + Xi) — p\. 
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The quantities to, x%, p\ and V\ are chosen such that f'{0) = and the function 

u (x) = f'(U (x)) = 4(2 - arctan(x + Xi)) 1/3 - v x (20) 

satisfies conditions @ corresponding to blow-up at t — x — with U — 0. Values of 
these quantities (with first three decimal digits) are 

to = -1.155, X! = 0.183, pi = 1.818, v x = 4.882. (21) 

Convergence ( T15"]) of the scaled wave profile U\(t,x) = Q\U{t,x) to the universal func- 
tion is demonstrated in Fig. [Tb. 

3 Renormalization group approach 

The universality just described can also be explained using the renormalization group 
approach. For this purpose, we consider Q\ in ((7j) as an operator acting in the space 
of solutions of the Burgers equation ([T|) with initial conditions satisfying @. Then 
the system dynamics can be seen as an action of Q\ combined with the scaling of 
time, space and state. This operator defines a differentiable group with the property 
G\iGx 2 — Q\i+x 2 - The universal function w(t, x) represents a stationary point (ITT]) of 
the renormalization group operator. Our analysis showed that this stationary point 
(more precisely, a set of stationary points w(t,x) determined by Q up to a scaling 
constant c) is asymptotically stable in the sense of Lyapunov (for A considered as 
"time"). From this property, the universality result (fTO]) follows. 

In the case of a general conservation law ffT2]) . the function U\(t,x) = Q\U(t,x) is 
a solution for a conservation law with a different flux function f\{U) = A 2 / 3 /(A~ 1 / 3 f/). 
Thus, Q\ is a renormalization group operator acting in the functional space (U, f) t— > 
(U\, f\) of solutions and fluxes. The universality ( IT5]) of the blow-up is explained by 
the fact that Q\ has the asymptotically stable stationary point (U, f) = (w(t, x), U 2 /2) 
corresponding to the universal solution ([9]) of the Burgers equation (JlJ). 

The role of the renormalization group operator G\ is similar to those in other, 
much more sophisticated theories. For the inner scale, x ~ ut ~ t 3 / 2 in ®, the wave 
dynamics is governed by the universal function, which is a stationary point of Q\. This 
is analogous, e.g., to the stationary point of the renormalization group operator, which 
determines critical phenomena in second-order phase transitions [UJ . At larger spatial 
scales, there is no universality and the solution depends on the initial condition u (x) as 
well as on the flux function f(u). This is analogous, in turn, to the phenomenological 
Landau theory of second-order phase transitions valid at larger (though still small) 
deviations of temperature from a critical value [9]. 

4 Blow-up as a solitary wave in wavenumber space 

Finite time blow-up implies rapid increase of solution in a short wavelength range. 
We will see now that the universality of the blow-up in x-space induces a solitary 
wave moving with constant speed to large wave numbers in logarithmic coordinates. 
Consider the Fourier transform in x-space of the solution, u(t, k) = J : x [u{t,x)\. Since 
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u(t,x) is real, we have u(t, —k) = u*(t, k), and we assume k > in the analysis that 
follows. For Fourier transformed functions, the relation (11 01) yields 

w(t,k) = lim u x (t,k), (22) 

A— >oo 

with the /c-space renormalization group operator 

u\(t, k) = g x u(t, k) = A 4/3 w(A~ 2/3 t, Xk). (23) 

Expression ( 123]) can be checked using Fourier transform of (J7|). Note that the function 
w(t, x) grows as w ~ x 1//3 for large x, see ([9]), and its Fourier transform (regularized by 
adding a small imaginary part to x) behaves as w ~ fc _4//3 for small k. 

Let us consider the wave profile u(t, k) transformed to logarithmic coordinates 

r = -log(*/*o), f = logfc, (24) 

which we denote by the same letter u(t,£). Here the blow-up at t = corresponds to 
r — > oo. Using (1231 . we find 



«a(t, = e 4a/3 u (t + ^a, £ + , a = log A. 



Because of (1221). we have 



(25) 



e 4a/3 w Qa, £ + a^j iu(0, as a oo, (26) 

where we put r = 0, and w(t, £) denotes the universal function w(t,k) written in 
coordinates (124T) . Now we substitute a by 3r/2 and write (1261) in the form 

u (r, £ + -r) -)• e~ 2T w(0,O as r ^ oo. (27) 



Therefore, u(t,^) forms a wave in the £-space for large r with the profile described 
by the universal function w(0,C,). This wave moves with the constant speed 3/2 and 
decays as e _2r . 

Fourier transforms of derivatives u^ n '(t, x) = d n u/dx n have the form u( n \t,k) = 
(ik) n u(t, k). Using logarithmic coordinates (r, £), derivations similar to (j24"|) - (l2T|) yield 

w (n) (t, £ + -> z n e (3n/2 - 2)T e n ^(0, as r -> oo. (28) 

Recall that w ~ /c~ 4/3 = e~ 4€/3 for f ->■ -oo (small fc). Thus, e n? w(0,£) ^ as 
£ — >■ — oo when ra > 4/3. In this case the wave derivative forms a solitary wave in 
the space (r, £), which moves with speed 3/2, grows exponentially as e ( 3n / 2_2 ) T , and has 
the universal shape determined by e n ^w(0,£). Formation of such a wave is shown in 
Fig. [2] for the second derivative, n = 2. Note that a similar wave describes a finite-time 
blow-up in shell models of turbulence, see [5]. 
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Figure 2: Profiles of the scaled second derivative e~ T \d 2 u/dx 2 \ in logarithmic Fourier 
coordinates (r, £) defined by (JUJ) for the solution shown in Fig. [D(a). Initial condition 
corresponds to r = 0, and the blow-up corresponds to r — > oo. The profile forms a 
solitary wave traveling with constant speed 3/2, which reflects the universality of the 
blow-up in x- and /c-spaces. The bold red curve shows the exact limiting shape of the 
wave determined by the universal function w(t,x). 

5 Conclusion 

We showed that solutions of conservation laws in one space dimension have universal 
structure when approaching the finite-time singularity (blow-up). The limiting wave 
profile is described by a scaled universal function in x-space, or by a solitary wave 
moving with constant speed in logarithmic coordinates of (wave number) fc-space. The 
universal function is independent both of initial conditions and of specific expressions 
for fluxes of conserved quantities. The universal properties are confirmed numerically 
in Figs. [U and El Experimental observation does not seem to be difficult due to large 
variety of nonlinear waves described by conservation laws in physical system. The 
obtained results can be extended to higher dimensional spaces, in particular 2D and 
3D, where the universality can be used as an effective tool for identification of blow-up 
singularities. This extension will be published elsewhere. 
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